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I. INTRODUCTION 

Coupled quantum dots (QDs), often referred to as 
quantum dot molecules (QDMs), have recently attracted 
much attention^ due to their potential application 
in various implementations of quantum computation 
schemes. Specifically, there have been many proposals 
for employing two-electron spin states in QDMs 3 -^ 5 -, ben- 
efiting from long decoherence times of the spin 6 . For 
instance, it was suggested to use singlet and triplet 
states as logical qubit states and to perform quantum 
computation 7 and Bell-states measurements 8 . Initial- 
ization, control and read-out of the state of two con- 
fined electrons in a QDM have already been experimen- 
tally demonstrated 2 -. Moreover, such structures are pro- 
posed for coherent optical manipulation of two-electron 
states^. 

Semiconductor QDMs are embedded in a solid state 
environment, which leads to electron interaction with 
the phonon reservoir. The presence of phonon-mediated 
coupling between energy eigenstates of a QDM leads 
to new effects in the physics of these structures, as 
compared to individual QD a 11 ' 12 . In particular, if the 
lowest states correspond to electron localization in two 
different dots, the relaxation between these states has 
the character of phonon-assisted tunneling, consisting in 
the transfer of an electron from one dot to the other. 
Such a process results from an interplay between the 
carrier-phonon coupling and tunneling coupling between 
the dots, which is a desirable element of many pro- 
posals of QDM-based quantum computing. Phonon- 
assisted tunneling has been thoroughly studied in the 
case of QDMs doped with a single electro n 13 ' 14 ! 15 ' 16 . Also 
phonon-induced triplet-singlet relaxation (via spin-orbit 
coupling) has been analyzed^ 7 -. However, to our knowl- 
edge, spin-conserving relaxation between singlet states of 
a two-electron system has not been discussed. 

In this paper, we analyze phonon-assisted tunneling 
in a quantum dot molecule consisting of two laterally 



coupled semiconductor quantum dots. A system doped 
with two electrons is considered. We study singlet-singlet 
relaxation, that is, relaxation between the two lowest 
states of two electrons in a QDM corresponding to the 
singlet spin configuration. For a specific GaAs QDM sys- 
tem, phonon-assisted relaxation rates are calculated. As 
we will show, in the parameter areas where the relax- 
ation is efficient, it involves charge transfer between the 
dots. Thus, it represents a phonon-assisted tunneling 
process. We study how the Coulomb interaction in the 
two-electron system influences the relaxation of electrons 
in comparison with the case of a QDM doped with a sin- 
gle electron. It is demonstrated that the presence of one 
electron strongly affects the tunneling of the other. As a 
result, the rates of the phonon-assisted electron tunnel- 
ing for the two doping cases (with one or two electrons) 
differ considerably, which is especially apparent in their 
dependence on the distance between the constituent QDs. 
We consider electrons interacting with acoustic phonon 
modes via deformation potential and piezoelectric cou- 
plings and show their distinguished impact on tunneling 
in QDMs. It is shown that the piezoelectric mechanism 
resulting from the considerable change of charge distri- 
bution is of great importance in the considered system 
and for some ranges of QDM parameters it is even the 
dominant contribution to relaxation. We show that the 
phonon-assisted tunneling is strong in coupled quantum 
dots and one should be aware of its influence when de- 
signing quantum computation schemes in QDMs. 



The paper is organized as follows. In Sec. [TTJ we intro- 
duce the model describing a quantum dot molecule with 
the Coulomb interaction and coupling to the phonon en- 
vironment. Section IlIII contains the results on phonon- 
assisted tunneling rates for the two systems under con- 
sideration. In Sec. IIV1 we conclude the paper with final 
remarks. In the Appendix, we summarize the theory of 
single-electron phonon-assisted tunneling. 



2 





-20 20 40 60 
x (nm) 



FIG. 1: (a) Schematic plot of a laterally coupled double quan- 
tum dot. (b) Lateral confinement potential describing the 
double QD structure. 



II. MODEL 

A. Electron states 

We consider a quantum dot molecule which consists 
of two laterally (in x direction) coupled quantum dots 
[see Fig. HJa)]. The structure doped with two electrons 
is studied. The Hamiltonian of the electron subsystem is 
given by 



H e = 



2m' 



(V 2 + Vg) + C/(r a ) + U(r h ) + V(v a , r b ), (1) 



where m* — 0.07m is the effective mass of an electron 
in GaAs. C/(r a / b ) is the confinement potential for two 
electrons referred to as 'a' and 'b', respectively, and 



V(r a ,r b ) 



1 



47re £ r 



is the Coulomb interaction between the electrons. Here, 
e denotes electron charge, Sq is the vacuum dielectric 
constant, and e r is the static relative dielectric constant. 
We assume a separable confinement potential 



U{r) = U(z) + U(y) + U(x) 



(2) 



\m*iol z 1 + \m*uly 2 + U(x), 



where U{z) is the potential in the growth direction and 
U(y) is the lateral confinement potential. The potential 
describing the double quantum dot structure is U{x) and 
has two minima, defining the two QDs. We choose it in 
the form: 



U(x) 



- Uq exp 



1 /x-d/2 



(3) 



-(U + AU) exp 



1 fx + d/2 
2 



This model potential has the advantage that it is smooth 
and allows one to independently control the distance be- 
tween the dots d and the depths of both potential wells, 



Uq and Uq + AU. The difference between the depths of 
the two constituent dots, AU, is referred to as the offset. 

The dynamics in the growth and lateral y directions 
is restricted to the respective ground states, which are 
described by Gaussian wave functions 



1 



v{y) = 



i 



exp 



z 

2h? 



(4) 



(5) 



Here, h denotes the electron wave function width in the 
growth direction z, while / is the width in the lateral 
direction y. The restriction to the ground states is a rea- 
sonable assumption in the considered confinement con- 
ditions, since the energy separation from the next eigen- 
states has typical values larger that 10 meV so that these 
states do not contribute to the studied dynamics. The 
complete wave function of a single electron can be written 
in a product form 

*n(r)=it> n (x)<p(y)<l>(z), (6) 

where ip n (x) is the nth lowest state of electron in a QDM 
obtained together with its eigenenergies from the numer- 
ical solution of Schrodinger equation. The two lowest 
single particle eigenstates are described by the wave func- 
tions $o( r ) and < I ) i(r) with the corresponding energies eo 
and t\. 

In order to analyze the relaxation mechanisms for 
a system doped with two electrons, we construct two- 
particle spin-singlet states 

|RL) = -j= ^aoT a U + a iT a oj.) l vac )> ( 7 ) 

|RR) = aJ T 4jvac). (8) 

Here, aj,^,^ creates an electron in the Oth (1st) lowest 
single-particle state with spin up (spin down), and |vac) 
denotes an empty quantum dot system. The correspond- 
ing spatially symmetric wave functions are 



*RL(r a ,rb) 
*RR(r a , rb) 



$o(r a )$i(r b ) + $i(r a )$ (r b ) 



V2 



$ (r a )$ (r b )- 



(9) 
(10) 



For the considered two-electron system, we include the 
Coulomb interaction between electrons and solve the sec- 
ular equation in the subspace spanned by the states |RL) 
and |RR), with the projected Hamiltonian 



H = 



eo + ei + vqo Vqi 



v i 2e Q + Vn 
where the Coulomb matrix elements are 



voo = V J d 3 k ^{Rc[^* (fe)^ 11 (fc)] + |^oi(fe)| 2 }, 
voi = V2V J d 3 k-^Re[T* (k)T i(k)], (11) 
vn = V J d 3 k ^|^ 00 (fc)| 2 , (12) 
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with 



Vo 



8n 3 eoe r a 

The single-particle form factors are defined as 



We consider the electron relaxation in the double-QD 
structure, which is a real transition on a picosecond time 
scale, therefore it can be treated within the Fermi golden 
rule approach. The coupling between the two considered 
singlet states in a two-electron configuration, resulting 
from the carrier-phonon interaction Hamiltonian (|15p is 



T nm {k) = J d 3 r $;(r) e ' ;fer <Mr) (13) H int = |0)(1| £ G s (k) (& >ifc + &+,_ fc ) +H.c. 



and for our choice of Gaussian wave functions in the y 
and z directions are 



T nm {k) = cxp z - — 



4 4 



(14) 



The resulting eigenstates of the interacting system are 
labeled as |0) and |1) and the corresponding wave func- 
tions can be written in the form 

a T • a T 

*o = cos-Wrl +sin-WRR, 

• a T a -r 

*i = -sm-VRL + cos-WRR, 



where 



a = arctan 



''01 



£o - ei + uix - woo 



and the energies are Eq and £i, respectively. The split- 
ting between the two-electron energies is 



AE = \E X - E \ = J (e - e x + «n - v oa ) 2 + Av 2 



oi- 



B. Carrier-phonon coupling 

In the considered QDM system, carriers not only in- 
teract with each other, but are also coupled to phonons. 
The free phonon Hamiltonian is 

where b e k and b s k denote phonon creation and annihila- 
tion operators, respectively. The corresponding frequen- 
cies are u) s ,k, where s labels different phonon branches 
and fe is the phonon wave vector. 

The interaction of the electrons with the phonon reser- 
voir is described by the Hamiltonian 



Hi 



in. = °« 

nm,er 



aa m ,„J2Fs,nm(k)(b s<k + bl_ k ), (15 



where F s ^ nn i (fe) are the single-particle coupling constants 
[see Eqs. (IA.1[) and (|A.2|) in the Appendix], which have 
the symmetry F snn t(k) — F* n , n (— fe), and a denotes the 
spin orientation. 



where G s (fe) are the two-electron coupling constants 
(given below). 

The energy difference between the electron states is 
considerably smaller than the energy of longitudinal op- 
tical phonons (LO), which is 36 meV in GaAs, thus they 
will not contribute to the relaxation mechanisms. There- 
fore, we consider interaction only with the relevant acous- 
tic phonons via the deformation potential and the piezo- 
electric coupling. 

Using the carrier-phonon coupling constant for the de- 
formation potential interaction [Eq. (|A.1[) ]. one finds the 
effective coupling between the two-electron states 



GP P (fc) = 



/ hk 

2pVc x 



D c g{k), 



where p is the crystal density, V is the normalization 
volume of the phonon modes, c\ is the longitudinal speed 
of sound, and D e is the deformation potential constant 
for electrons. The geometrical properties of the wave 
functions are reflected in the form factor G(k), which has 
the form 

G(k) = J d^ J d 3 r h %(r ai r h ) (16) 

x (e lfcr »+e lferb )*i(r a ,r b ). 

It can be written by means of the single-particle form 
factors defined by Eq. (fl3|) as 

Q(k) = ^sina^ooCfe) -^ii(fe)] +V / 2cosa^ r oi(fe) 



( k 2 h 2 k 2 l 2 \ - 
exp --£- S_ g {kx ). 



(17) 



The two-particle coupling elements for the piezoelectric 
interaction arc 



G* E (fe) = -u 



dpe 



2pVc s k £ £r 



M s {k)g{k), (is) 



where c s is the speed of sound (longitudinal c\ or trans- 
verse Ct , depending on the phonon branch) and dp is the 
piezoelectric constant. The function M s (k) does not de- 
pend on the value of the phonon wave vector, but only 
on its orientation. For a zinc-blende structure, it reads 



Ms{k) — k x (e Si k)yk z + [s St k)zky 



(e^.fc) zk x + {Cs,k)xk 2 



~\~k z {c s k ) x ky -(- i^c s k )yk x 



(19) 
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where e St k is the unit polarization vector for the phonon 
wave vector k and polarization s, and fc = fe/fc. We 
choose the following phonon polarization vectors 



e\k — k = (cos sin# cos <p, sin# simp) , 
eti.fc = (0, sin </>,- cos (p) , 
&t2,fe = ( — sin cos $ cos (p, cos sin ip~) , 

for which the functions M s (fe) read 
3 



(20) 



Mi(fc) 



• sin#sin(2#) sm(2<^), 



(21) 



M t i(fe) = -sin(26>)cos(2^), 
M t2 (fe) = sin0(3cos 2 6>-l)sin(2^). 

The properties of the phonon environment are repre- 
sented by phonon spectral density 



R(u) = ^|nH + l|^|G s (fc)| 2 



x [(5(w - w s ,fc) + <S(w + w s ,fe)] , 



(22) 



where n(w) is the Bose distribution function. The defor- 
mation potential contribution is 



i? DP H 



i^V |n(w) + l| 

I dip sinO d9 \g(u> cos 6/ ci)\ 2 



(23) 



x exp 



— j sin 2 9 {l 2 cos 2 ip + h 2 sin 2 tp) 
2ci 



where 



nDP _ 

H - 



I? 2 



The piezoelectric term is 



167T 3 ft/9cf 



2 k- 



i(> |n(w) + l| / d<£ / sin 9 d9 
Jo Jo 

x\M s (k( v ,9))\ 2 \G(uj cos 9 /c s )\ 2 

2 

2 /) A2 „„„2 , j 2 • 2 



x exp 



—j sin 2 6> (Z 2 cos 2 + ft, 2 sin 2 tp) 
zee 



where 



pPE _ 



167T 3 /ipC 3 £g£ 2 



Note that the coupling constants for deformation poten- 
tial and piezoelectric channels have different parity (as 
functions of fe) so that these two transition channels do 
not interfere. 

In order to study phonon-assisted relaxation, we em- 
ploy the Fermi golden rule and obtain the relaxation rate 



/AE\ 



(24) 
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i 


10.0 nm 


s-direction 
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TABLE I: The GaAs material parameters and QDM system 
parameters. 



which is proportional to the phonon spectral density at 
the frequency corresponding to the splitting energy AE. 

The material parameters (corresponding to GaAs 
quantum dots) and parameters of the QDM system are 
given in Table [U Moreover, details concerning relaxation 
in a single electron QDM system, which will be used for 
comparison, are presented in the Appendix. 



III. RESULTS: TUNNELING RATES 

In this section, the results for phonon-assisted transi- 
tions in a double quantum dot are presented. We con- 
sider a QDM doped with two electrons and the singlet- 
singlet relaxation channel. In order to investigate how 
the Coulomb interaction influences the relaxation pro- 
cesses, the results are compared with those for a sin- 
gle electron case, calculated in a way similar to Refs. 
Il3lfl4lfl5lfl6l (see Appendix). The quantitative results are 
obtained at temperature T = K for GaAs quantum 
dots with the sizes h = 4 nm and I = a = 10 nm in 
growth and lateral directions, respectively. 

The probability of phonon-assisted electron transitions 
[Eq. (|24p] is proportional to the spectral density of the 
phonon reservoir at the frequency corresponding to the 
energy splitting AE. Therefore, the transition rate will 
be high when this energy lies in the frequency range of 
maximal values of the phonon spectral density. In order 
to see which parameter range is favorable for relaxation, 
we first study the energy splittings and phonon spectral 
densities for the two considered doping cases. 

In Fig. [2a), we plotted the energy splitting for a sin- 
gle electron in a QDM as a function of the confinement 
depth offset AU for Uq = 30 meV and a few values of the 
distance d between the QDs. The minimum value occurs 
always when the QDs are the same and, in general, is 
smaller for larger distances between dots, where they do 
not influence each other. In the case of two electrons in 
a QDM [Fig. Hfb)], the splitting energies have a slightly 
more complicated behavior. Now, the minimum value 
is shifted due to interplay between the on-site (single- 
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FIG. 2: (a) Splitting energy as a function of the QDs offset for 
a single electron in a QDM for different distances d between 
QDs. (b) As in (a) but for two electrons in a QDM. 
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FIG. 3: (a) Total spectral density of the phonon reservoir for 
one electron in a QDM for AU = meV and different dis- 
tances d between the QDs. (b) As in (a) but for two electrons 
and AU = 11 meV. 

particle) potential and the Coulomb interaction, which 
also depends on the distance between QDs. The split- 
ting energies are larger, since they describe two-particle 
states affected by the Coulomb coupling. While for one 
electron, the energies are symmetric with respect to the 
resonance point (minimum splitting), in the two-electron 
case this symmetry is lost, except for very large separa- 
tions. 

Since the wave functions obviously depend on the con- 
finement offset AJ7, so do the coupling constants and, in 
consequence, the phonon spectral densities. In order to 
gain some information on their general shape, we study 
the spectral densities for the values of the offset which 
correspond to the minimal values of energy splitting. In 
Fig. [U(a) , the phonon spectral density for a single elec- 
tron [see Eqs. (|A.3[) and (|A.4|1 in the Appendix] is plot- 
ted for the offset AU = meV. In general, the values 
of phonon spectral densities depend on the overlap be- 
tween the wave functions and thus are large for small 
distances d between the QDs. The function has its max- 
imum for lo w 0.4 meV and a cut-off at u) rs 2.5 meV. 
One can expect high rates for energy splittings from 0.1 
to 1 meV, especially for small distances d. From Fig. [U(a) 
it is clear that for closely spaced QDs, the energy splitting 
is larger than 2 meV and lies almost beyond the cut-off 
of the phonon density, which will result in lower transi- 



FIG. 4: (a) Total phonon spectral density with the contri- 
butions resulting from deformation potential and piezoelec- 
tric couplings for one electron in a QDM for AU = meV 
and d = 38 nm. (b) As in (a) but for two electrons and 
AU = 11 meV. 



tion rates. For larger distances the splitting is smaller, 
but also the amplitude of the spectral density is smaller. 
The interplay of phonon density and splitting energies 
will be reflected in nontrivial dependence of relaxation 
rates on the distance between the QDs. 

For a two-electron QDM, the phonon spectral density 
has, in general, smaller values [Fig. [3jb)], since the over- 
lap between corresponding two-electron wave functions 
is smaller. In this case, the cut-off energy (ui w 3 meV) 
as well as the energy splitting is larger. One can see that 
phonon- assisted transitions in both systems will be large 
for energy splittings smaller than 3 meV and will strongly 
depend on the distance d. 

The electron-phonon interaction via both deformation 
potential as well as piezoelectric coupling is considered 
next. In order to see which interaction has a stronger in- 
fluence, in Figs.BJa) and[4|b) we present the total spec- 
tral density of the phonon environment together with the 
two contributions for a fixed distance d = 38 nm. It is 
clear that piezoelectric coupling in double quantum dot 
structures is of great importance in contrast to optical 
processes in single QD structures, where this interaction 
can in many cases be neglecte d 18 ' 19 ! 20 . This results from 
the fact that electron relaxation induces a large change 
of charge redistribution, especially when it involves tun- 
neling to the other dot. Since in a single electron as well 
in a two-electron system the two phonon contributions 
may cover different frequency sectors, they will also play 
a role in the transition rates in distinct parameter areas. 

We start the discussion of phonon-assisted relaxation 
from the deformation potential contribution. For a one- 
electron QDM [Fig. Ufa)] , the rates are symmetric with 
respect to the offset of the quantum dots AU. This re- 
sults from the symmetric behavior of the splitting en- 
ergies AE. When the QDs are close to each other, 
d < 28 nm, the rate is low, since the energy difference be- 
tween the two lowest electron states is much larger than 
the acoustic phonon energies. Thus one-phonon transi- 
tions are impossible. For distances d from 28 to 38 nm, 
the relaxation rate is high. It corresponds to the situation 
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FIG. 5: (a) Electron relaxation rate assisted by phonons via 
deformation potential coupling for a single electron in a QDM 
as a function of QDs offset and distance d. (b) As in (a) but 
for two electrons. 



FIG. 6: (a) Electron relaxation rate assisted by phonons via 
piezoelectric coupling for a single electron in a QDM as a 
function of QD confinement offset AU and distance d. (b) As 
in (a) but for two electrons. 



when the energy splittings are comparable to the phonon 
energies. The transition rate reaches its maximum value 
of 30 ns^ 1 for d « 32 nm. Here, the relaxation conditions 
are most favorable, since the distance between the QDs is 
large enough for the splitting energy to coincide with the 
maximum value of the phonon spectral density. For large 
distances, d > 38 nm, the rate vanishes in spite of small 
splitting energies, since the overlap between the electron 
wave functions tends to zero and, in consequence, the 
spectral density vanishes. The transition rates are also 
small for large offsets, \AU\ > 3 meV, since it leads to 
large energy gap between the levels. 

For a two-electron QDM [Fig. [5^b)], the maximum of 
the relaxation rate shifts with growing distance towards 
larger confinement offsets, which was already visible in 
the splitting energies. Larger distances d between the 
QDs are needed for efficient relaxation, which is an evi- 
dence of the Coulomb interaction between two electrons, 
leading to an increase of the splitting energies. In general, 
the maximum magnitude of the relaxation rates is com- 
parable with that for a single electron but the parameter 
range in which their values are maximal is shifted due to 
the electron-electron interaction. 

In the case of piezoelectric coupling [Figs. [Ha) 
and EJb)] , the relaxation rate has relatively large values 
in a smaller range of QD offsets. This is a result of a dif- 
ferent form of corresponding spectral density, which, in 
general, is narrower than for the deformation potential. 
Therefore, smaller splitting energies are more favorable. 
For the same reason, it is shifted towards larger distances 
d. The relaxation rate reaches the values of 30 ns -1 , 
which is as large as that for deformation potential. This 
maximum appears at the distance d w 36 nm for a single 
electron and d w 39 nm for a double electron QDM. 

The total phonon-induced relaxation rate, including 
both deformation potential as well as piezoelectric con- 
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FIG. 7: (a) Relaxation rate for all phonon modes and one 
electron in a QDM. (b) As in (a) but for two electrons. 



tributions, is shown in Figs. EH a) and [7jb) as a func- 
tion of QD offset AU and separation d. For a sin- 
gle electron system, the rate is high for offsets between 
AU = —3 and AU = 3 meV and for distances from 
d = 28 to d = 42 nm, and reaches its maximal value 
of 35 ns -1 for identical QDs separated by the distance 
of d w 34 nm. In case of the QDM doped with two 
electrons, the relaxation mechanism is strong for offsets 
between AU ss 9 and AU « 13 meV and distances from 
d « 30 to d « 44 nm. Its maximum value also reaches 
35 ns -1 for AU « 11.5 meV and d « 36 nm. 

To understand the nature of the relaxation process in 
the two-electron case, in Fig. [8j we plot the difference 
between the average number of electrons in the left quan- 
tum dot in the two-electron singlet states. One can see 
that in the area of efficient relaxation (cf . Fig. [7]) , the av- 
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FIG. 9: (a) Total phonon-assisted relaxation rate with two 
contributions in a two-electron QDM for d = 35 nm. (b) As 
in (a) but for d = 38 nm. 



FIG. 8: The difference between the average number of elec- 
trons in the left quantum dot in the two-electron singlet 
states. 



IV. CONCLUSION 



erage electron number changes in most cases almost by 
one. This shows that the relaxation in the two-electron 
case is associated with a considerable charge transfer and, 
therefore, can be interpreted as a phonon-assisted tunnel- 
ing process. 

For a single electron, the energy eigenstates follow a 
universal model of level anticrossing, with the energy 
splitting AE = ^{AU) 2 + At 2 , where t is the "tunnel- 
ing matrix element" , corresponding to half of the mini- 
mum energy splitting in Fig. [2^. This element affects the 
phonon-assisted tunneling rate in a twofold way. First, 
it determines the splitting of the energy levels and its 
position with respect to the area of large phonon spec- 
tral density. Second, it affects the degree of mixing of 
the wave functions, thus directly changing the spectral 
density. It should be noted, however, that the relaxation 
rate cannot be fully characterized by this single param- 
eter, since the phonon spectral density depends on the 
actual geometry of the system, and therefore the spa- 
tial separation between the dots is itself of direct impor- 
tance. In the two-electron system, the situation is even 
more complicated, since the energies and wave functions 
are affected by the interplay between the single-particle 
"tunnel coupling" and the Coulomb interaction between 
the electrons. This is manifested in the increased reso- 
nance width and loss of symmetry in Fig. [2)3. As a result, 
the relaxation rates are also asymmetric with respect to 
AU [see Figs. [7(a) and E£b)]. 

To have a better insight into particular phonon con- 
tributions, in Figs. [5(a) and [DJb) we present the total 
relaxation rate along with the contributions from both 
the coupling mechanisms for a two-electron QDM for 
c? = 35 nm and d = 38 nm, respectively. For relaxation 
rates of comparable values, the dominant phonon cou- 
pling can be different. For instance, for d = 35 nm, the 
deformation potential coupling is crucial and the piezo- 
electric effect is a few times smaller, while for d = 38 nm 
the situation is reverse. 



In the present paper, we have studied phonon-assisted 
relaxation and tunneling in a quantum dot molecule. 
Structures doped with two electrons have been consid- 
ered and compared with the case of a single electron. By 
comparison of these two systems, it was shown that the 
Coulomb interaction influences the tunneling rates and 
leads to energy renormalization and shift of the range of 
efficient relaxation. We studied in detail carrier-phonon 
interactions via both deformation potential and piezo- 
electric coupling and showed the difference in their be- 
havior and impact on relaxation. We have shown that 
the relaxation in the two-electron case is accompanied 
by a charge transfer between the dots and, therefore, can 
be regarded as a phonon-assisted tunneling process. 

It should be noted that the values of phonon-assisted 
tunneling rates in a QDM system are comparable with re- 
laxation times in a single QD^i. Moreover, in comparison 
with the spin coherence times being up to milliseconds 6 , 
the phonon-assisted relaxation times are up to several or- 
ders of magnitude faster. This shows that while design- 
ing the quantum computer implementations on electron 
states in double quantum dots, one has to take into ac- 
count the coupling of the carriers to the phonon degrees 
of freedom. Finally, it should be noted that the calcula- 
tions were performed for zero temperature, which gives 
a lower bound for tunneling rates. 
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APPENDIX: SINGLE ELECTRON TUNNELING 

In this Appendix, the formalism for the tunneling in 
a single electron QDM system is presented. In this case, 
the electron Hamiltonian [Eq. (JXJ) ] is reduced to a simpler 
form, 



h 2 „ 

He = ~ V 2 

2m* 



and the electrons are described by a single particle wave 
functions given by Eq. ^ . We label the two lowest single 
electron states as |0) and |1). 

The relevant part of the carrier-phonon interaction 
Hamiltonian describing electron transitions between the 
constituent QDs is 

Hint = |OXl|X>.,oi(fc) (b., h + bl_ k ) +H.c, 



where the single-particle coupling constant for the defor- 
mation potential is 



' hk 

2pVc A 



D e T 01 (k), 



(A.l) 



with the form factors given by Eq. (|13p . The coupling 
element for piezoelectric interactions reads 



M,(k)F 01 (k), (A.2) 



where the functions M s (k) are given by Eq. (f2T|) . 

The corresponding phonon spectral densities for a sin- 
gle electron in a QDM are 



CM = iz£"V>M + i| 



(A.3) 



x / dip sin 9 dO \T i (locos 6/ ci)\* 



x exp 



■ sin 2 6 (l 2 cos 2 ip + h 2 sin 2 tp) 



and 



^ie sM = iZo* w |n(w) + 1| / dip / sinOde (A.4) 

Jo Ja 

x\M s (k(<p, 9))\ 2 \T i(iJcosd/c s )\ 2 



x exp 



to 



— sin^ 9 (r cos z <p + hr sin <p) 



where 



T {)l {k x ) = / dx r {x)e^ x ^{x) 



The Fermi golden rule relaxation rate is then calcu- 
lated from Eq. (j2"4"| . using the total spectral density in- 
cluding both relaxation channels. 
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